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The spin-down power of an isolated neutron star can drive its central density in- 
crease and overall structural changes, and trigger a quark-hadron phase transition. 
A series of observational signatures may be seen as a result of the phase transi- 
tion, including pulsar spin-down and glitch behaviors, Soft Gamma Repeaters or 
Gamma-Ray Bursts. 



1 Introduction 

The sky has long been a successful test ground and discovery site for fun- 
damental physics. Examples include the discovery of Helium frouuthe solar 
spectrum in 1868, evidence for gravity waves from binary pulsars u in 1978, 
and more recently in 1985, the discovery of Bucky Balls {G&q), though not 
directly in the sky, through experiments designed to simulate environments 
near stars in favor of forming long chain carbon molecules a. 

Quark-hadron phase transitions are being pursued in accelerators in the 
high temperature regime. Meanwhile, the densities inside neutron stars may 
be high enough to allow the existence of deconfined quark matter. 

1.1 A Neutron Star with a Quark Core 

The structure of a neutron star with a quark core can be solved using standard 
equations of stellar structure B together with the equations of state (EOSs) of 
quark matter and normal neutron matter. The EOS of quark matter is much 
softer than that of neutron matter because of QCD asymptotic freedom. A 
neutron star containing a quark core is thus more compact than a normal 
neutron star. Consequently, the star has a larger maximum spin frequency. 
However, submillisecond pulsars have not been discovered and hence quark 
stars have not been identified this way. 

1.2 Rotating Relativistic Stars 

To solve the structure of a rotating relativistic star, Hartle developed §@ a 
perturbation solution based on the Schwarzschild metric of a static, spherically 
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symmetric object. Rotation distorts the star away from spherical symmetry, 
and its perturbed metric has the form 

ds 2 = ~e 2v dt 2 + e 2X dr 2 + e 2 ^d6 2 + e 2 ^(d</> - udt) 2 + 0{n 3 ). (1) 

In this line element, uj is the angular frequency of the star's fluid in a local 
inertial frame and depends on the radial coordinate r. It is related to f2, 
the rotational frequency seen by a distant observer. In this w.ay,|-the mass- 
central density relation can be solved for rotating compact stars &B0. We note 
that, the mass "increase", which is a main result in these works, represents 
the mass difference between two stars at different angular velocities with the 
same central density. Hence, these theoretical results are about star families 
relating to observations through pulsar statistics, and are difficult to be tested. 

2 Predictions on Quark-Core Neutron Star Spin-Down 
Behavior 

2.1 Another Way to Apply Hartle's Perturbation Method 

Instead of deriving mass-central density relations for a set of neutron stars, 
we have suggested to trace the evolution of an isolated neutron star during 
its spin-down process p. Associated with the central density increase, neu- 
tron matter is continuously converted to quark matter. The overall structure 
and spin-down behaviors of the star are modified. In Hartle's perturbation 
method, the central density is fixed as an input parameter, then stellar masses 
are solved as a function of angular frequency. In principle, we can derive a 
set of equations parallel to those of Hartle, with a fixed mass as an input and 
central density as a perturbation, and solve the central density and stellar 
structure as a function of angular velocities. In practice, it is more convenient 
to apply Hartle's method in a different way to solve the change of central 
density of an isolated star. 

An isolated neutron star has approximately constant gravitational mass 
during its spin-down process. The binding energy of a typical neutron star 
ly^lO 53 erg, and the ratio of rotational energy to gravitational energy is less 
than 0.1 for the most rapidly rotating neutron star 0. Since M Q ~10 54 erg, 
the rotational energy is only 1% of the stellar mass. In calculating the stellar 
configurations at different angular velocities, we need first to plot the mass- 
central density relatfon at different angular velocities using Hartle's method 
as in previous work aau. Under our approximation, we cross these curves with 
a line of constant mass, and the corresponding central densities are those of 
a rotating star at different angular velocities. 
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2.2 Pulsar Glitches 



The quarks inside the star are charged particles and are coupled with the 
solid crust through magnetic fields. They rotate at a different velocity from 
superfluid neutrons. During the spin-down process, the central density and 
the quark core grow and so does the fractional moment of inertia of quark 
matter as compared to the whole star (I q /I). In a simple starquake model, 
the post-glitch behavior of a pulsar can be described approximately as u, 

n(i) = f7 (t) + Af7 [Qe-*/ T + l-Q], (2) 

where Q is the healing parameter of the pulsar. Q = — TAfl(t=0)/ Aft, and 
r = — Af2/Af2(i=0). It has been shown that QocI n /I (ref. 3), where I = I n + 
I c is the total moment of inertia of the star, /„ is the moment of inertia carried 
by neutral particles, and I c is that of charged particles. Hence, a decreasing 
Q over a long term may indicate continuous conversion from neutron matter 
to quark matter inside the pulsar. 

2.3 Spin- Down History 

The pulsar spin-down can be described as i 

n = -m n , (3) 

and the braking index n increases as the moment of inertia decreases, 

An/n oc -AI/I. (4) 

As a result, |fi| is larger because the quark-core neutron star has a smaller / 
and thus is easier to brake. At the same time, the star tends to "spin-up" 
during the early stages of the phase transitions causing a higher £1 than in the 
normal neutron star case. Consequently, the "characteristic age" 

t = n/2\h\ = k- 1 n 1 - n (5) 

may underestimate the true age of the neutron star with a quark core i, while 
for a normal neutron star r is an upper limit. 

3 Onset of the Phase Transition: Catastrophic or Continuous? 

3.1 Event Rate 

In the previous section we considered the neutron star with a developed quark 
core in it. What if the central density of a normal neutron star happens to 
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increase from sub-critical density to critical density? Apparentljz-jbhe chance 
for this to happen is low, and we can estimate the event rate as Ii3 

where -Rns is the average neutron star birth rate in units of per year per galaxy 
and Pi is the initial spin period. R does not strongly depend on the exact 
critical density. The event rate is low because the central density increase is 
small if the initial spin is slow. However, equation @ suggest that these types 
of events happen at a finite rate, especially when we have millions of galaxies in 
our view. While the accretion power in astrophysics has been stressed in many 
astronomical environments, equation (||) suggest that the "spin-down power" 
is not negligible. The spin-down power can apply to any rotating systems with 
a critical point, not limited to quark- hadron phase transitions, but include 
Kion or Pion condensations inside neutron stars. Even for a rotating white 
dwarf the spin-down power alone may trigger a type la supernova as the 
Chandraseckhar mass limit depends on the rotational velocity. 



3.2 Catastrophic 

If the phase transition is catastrophic, it may happen on a time scale of sec- 
onds. In such a short time the star contracts from a radius of ~15 km to ^10 
km, associated with a sudden spin up. Most importantly, the gravitational 
energy released 

_ GM 2 fAR\ in53 fAR\ 
may be large enough to power a cosmological gamma-ray burst (GRB) 



3.3 N on- Catastrophic 

Even though quark-hadron phase transition is likely to be first order, it is dif- 
ferent from water-vapor phase transition in that the former has an additional 
freedom of whether a proton deconfines to uud quarks or a neutron becomes 
udal quarks, which does not conserve electric charges in the two phases al- 
though the overall charge neutrality is conserved with the help from leptons. 
Hence, the two phases are not necessarily separated by gravity EJ. In this 
case, the phase transition may happen slowly on a time scale of 10 5 years and 
the energy is released at an average rate of 10 40 ergs s _1 . The majority of 
this energy is released via neutrino emission, and only a tiny fraction is used 
to heat the star up to a surface temperature of 3x 10 6 K, yielding a soft X-ray 
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luminosity of ^10 35 ergs s _1 . Note that it is 25 tiroes more luminous than 
the sun, and it is not necessarily in a binary system 113. While the fluid core 
of the star is contracting, the solid crust may have stress built up in it, and 
the cracking of the crust can release bursts of energy 13. 

4 Observations 

In this section we confront theoretical predictions with observations. Radio 
pulsars constitute only a small fraction of all neutron stars because of the 
beaming nature of their radio emissions, and hence are not the best place 
to look for the signatures of the phase transition. So far only 3 pulsars have 
reliably measured n. Five glitches for Crab pulsar, and 7 for Vela are recorded 
with measured Q. Many more are needed to see the continuous phase transi- 
tion in radio pulsars. 

4-1 Soft Gamma Repeaters 

Soft Gamma Repeaters (SGRs) are X-ray transient sources (with ^20 keV 
photons) associated with young (10 4 yr) supernova remnants (SNRs). They 
are usually also quiescent X-ray emitters (with fcT~l keV, Lx^lO 35 ergs s" 1 ). 
So far 4 SGRs have been discovered in the Galaxy, and 1 in the Large Mag- 
ellanic Cloud. Two SGRs show fi/2|f2| (~10 3 yrs) < SNR age (~10 4 yrs) U. 
Since many more Anomalous X-ray Pulsars (AXPs) share some properties of 
SGRs, it seems that our prediction in previous section that wc should be able 
to see one at a time in our Galaxy is an underestimate. Note however that if 
the initial spin of a neutron star is 2 ms rather than 20 ms (that of the Crab 
pulsar), we get 100 times more event rate of phase transitiorL_Fast initial spin 
has indeed been suggested to produce strong magnetic fields t 3 ] and slow down 
neutron stars much faster than Crab- type pulsars. Radio astronomers may 
have been biased toward Crab-type pulsars when picking up neutron stars. 

4-2 Gamma-Ray Bursts 

GRBs are energetic explosions with total energy 10 51 — 10 54 erg (if isotropic), 
most of which is released with 100 keV to a few MeV photons. The time 
scales range from milliseconds to minutes. The event rate is estimated to be 
10 -5 yr -1 galaxy" 1 if the emission is isotropic, and can be much higher if the 
emission is highly beamed. We note that the event rate in equation (^) can 
still account for the GRBs even if they are beamed, and the faster initial spins 
of the neutron stars naturally offer large angular momentum for the sake of 
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beaming the electromagnetic radiations. 

Gravity wave detectors such as LIGO may be able to observe gravity waves 
associated GRBs, and tell which is the correct model. However, the gravity 
wave is emitted more or less isotropically while the electromagnetic radiation 
is likely to be highly beamed (with a beaming factor 100-1000) in the angular 
momentum direction. Hence, LIGO may need to accumulate 100-1000 events 
before seeing a gravity wave signal at the same time of a GRB. 

5 Conclusions 

The accretion power in astrophysics has been stressed over decades, while the 
"spin-down power" is often neglected. Here we show that, in any astrophysical 
systems with angular momentum and a critical point, the spin-down power 
can drive the system to the critical point at a small but finite probability. 
If the transition to critical state is energetic enough to be seen from distant 
galaxies, we can observe these events regularly, as in the case for the signatures 
of quark deconfinement phase transitions in radio pulsars, SGRs, and GRBs. 

Acknowledgments 

I thank Peter Williams for help with the manuscript. 
References 

1. J. H. Taylor, L. A. Fowler and P. M. McCulloch, Nature 277, 437 (1979) 

2. H. W. Kroto et al, Nature 318, 162 (1985) 

3. S. L. Shapiro and S. A. Teukolsky, Black Holes, White Dwarfs, and Neu- 
tron Stars, (Wiley, New York, 1983) 

4. J. B. Hartle, Astrophys. J. 150, 1005 (1967) 

5. J. B. Hartle and K. S. Thorne, Astrophys. J. 153, 807 (1968) 

6. F. Weber and N. K. Glendenning, Astrophys. J. 390, 541 (1992) 

7. M. Colpi and J. C. Miller, Astrophys. J. 388, 513 (1992) 

8. F. Ma and S. Luo, Bull. American Astron. Soc. 28, 961 (1996). 

9. F. Ma, Master's thesis, (Beijing Normal University, 1994) 

10. F. Ma, & B. Xie, Astrophys. J. Lett. 462, L63 (1996) 

11. N. K. Glendenning, Phys. Rev. D 46, 1274 (1992) 

12. F. Ma, Phys. Rev. Lett. 81, 1983 (1998) 

13. R. C. Duncan and C. Thompson, Astrophys. J. Lett. 392, L9 (1992) 

14. C. Kouveliotou et al, Astrophys. J. Lett. 510, LI 15 (1999) 



wsllwi: submitted to World Scientific on February 1, 2008 6 



